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ABSTRACT 

We present 35 ks Chandra ACIS observations of the 42 Myr old radio pulsar PSR B1451— 68. A 
point source is detected 0.32"±0.73" from the expected radio pulsar position. It has ^ 200 counts 
in the 0.3-8 keV energy range. We identify this point source as the X-ray counterpart of the radio 
pulsar. PSR B1451— 68 is located close to a 2MASS point source, for which we derive 7% as the 
upper limit on the flux contribution to the measured pulsar X-ray flux. The pulsar spectrum can 
be described by either a power- law model with photon index F = 2.4^03 and a unrealistically 
high absorbing column density iVn = (2.5^^ 3) x lO^-'^cm"^, or by a combination of a kT = 
0.35t!^:j7 keV blackbody and a P = 1.4 ± 0.5 power-law component for N^^ = 2.6 x lO^" cm^^ 
estimated from the pulsar dispersion measure. At the parallactic, Lutz-Kelker bias corrected 
distance of 480 pc, the non-thermal X-ray luminosities in the 0.3-8 keV energy band are either 
■^oT-skcV = (11-3±1.7) X lO^^ergs-i or LS°3-8keV = i^-^-tl) x lO^^ergs"!, respectively. This 
corresponds to non-thermal X-ray efficiencies of either '7o "a^^iskeV = ^oT-skcv/^ ~ 5 x lO^'^ or 
3 X 10^'^, respectively. 

Subject headings: X-rays: stars. Stars: neutron, pulsars: individual (PSR B1451— 68) 
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1. Introduction 

Old (> 1 Myr), rotational-powered pulsars have 
lost a siginificant amount of their initial rotation 
energy. Their spin-down power, E < 10"^^ erg s^^, 
is usually several orders of magntiude lower than 
those of younger pulsars. The lower energy budget 
translates into fainter high-energy emission. Only 
the closest of the old rotational-powered pulsars 
can be detected in X-rays. Since the bulk of the 
neutron star surface is too cold to be visible in 
X-rays (Yakovlev & Pethick 2004), one expects to 
see mainly the X-rays from the non-thermal mag- 
netospheric emission with a possible small thermal 
contribution from polar caps heated by infalling, 
accelerated particles (Harding & Muslimov 2001, 
2002). So far, eight old pulsars have been detected 
at X-ray wavelengths (e.g., Becker et al. 2004; Za- 
vlin & Pavlov 2004; Tepedelenlioglu & Ogelman 
2005; Zhang et al. 2005; Becker et al. 2006; Kar- 
gahscv et al. 2006; Misanovic et al. 2008; Hui & 
Becker 2008; Pavlov et al. 2009), several of them 
having only a few dozen counts. Most of the X-ray 
spectra can be described by relatively soft power 
laws with photon indices 2 < F < 4. Some of 
these spectra can be fitted with blackbody (BB) 
models or by a combination of BB and power law 
(PL) components. The BB models usually indi- 
cate a rather small emission area compared to the 
conventional polar cap sizes whose radii are cal- 
culated as Rpc = [27ri?^g/(cP)]^/^, assuming a 
magnetic dipole field (Harding & Muslimov 2001, 
2002). Here, i?NS is the radius of the neutron 
star, P is the pulsar period and c is the speed of 
light. There has been a debate as to how signifi- 
cant the thermal contributions arc (e.g., Becker et 
al. 2004; Zavhn & Pavlov 2004; Misanovic et al. 
2008; Hui & Becker 2008). Better count statistics, 
the consideration of new information from other 
wavelengths as well as a larger object sample are 
required to understand the quantitative contribu- 
tions of magnetospheric and thermal components. 

PSR B1451— 68 is an old, nearby pulsar, pre- 
viously unexplored in X-rays. The 0.263 s radio 
pulsar has a characteristic spin-down age r = 
P(2P)^^ — 4.25 X 10'' yr, a dipole magnetic field 
strength of 1.63 x 10^^ G, and a spin-down power 
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E = 2.1x 10^^ ergs ^. The proper motion compo- 
nents of PSR B1451— 68 were measured by Bailes 
et al. (1990): /i„ = -39.5 ± O.Tmasyr"! and 
/i5 = —12.3 ± 0.6masyr~'^. Unlike most other 
pulsars, the proper motion is directed toward the 
Galactic plane. Bailes et al. (1990) measured an 
annual parallax of tt 2.2 ± 0.3 mas and reported 
a distance of d = 450 ± 60 pc. According to Ver- 
biest et al. (2010) the distance corrected for the 
Lutz-Kelker bias is d = 480^60 pc- The spin-down 
power of PSR B1451— 68 is comparable to those 
of other nearby old pulsars which have been de- 
tected in X-rays. Aiming for constraints on the 
X-ray spectral properties of PSR B 145 1—68 , we 
report here on Chandra observations of this pul- 
sar. 

2. Observations and Data reduction 

PSR B 145 1—68 was observed with the Chan- 
dra X-ray observatory on 2010 May 16 using the 
ACIS-S detector in the VFAINT imaging mode. 
The nominal exposure time was 35.1 ks. The 
data reduction was done using CIAO (version 4.3). 
We reprocessed the data in order to apply the 
VFAINT background cleaning. We investigated 
event files with and without event position ran- 
domization using the Energy-Dependent Subpixel 
Event Repositioning (EDSER) algorithm by Li et 
al. (2004), as well as the subpixel algorithm by 
Mori et al. (2001). If not stated otherwise, we 
will discuss in the following the events with the 
pixel position randomization turned off. No large 
background flares occured during the observation, 
therefore we filtered the reprocessed event file for 
good time intervals (GTIs) which have less the 5<j 
deviation from the overall lightcurve mean rate. 
After the GTI filtering the actual exposure time 
was 34.7ks. Due to the ACIS frame time of 3.2s, 
it is not possible to search for X-ray pulsations 
corresponding to the P = 0.263 s of PSR B1451- 
68. The wavdetect task was applied to obtain 
an X-ray source list. A source was detected at 
the pulsar position on the back-illuminated chip 
S3. Using the CIAO task arfcorr we estimate 
that 95% of the source counts are included in a 
r = 2" aperture in an energy range from 0.3 keV 
to 5 keV. The aperture-corrected source count rate 
is 0.0059 ± 0.0004 cps. 

Since there is a stellar neighbour with small 
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Table 1: Chandra source positions in the repro- 
jected event files. 

method R.A. Dec. 

[°] ^ ^ [°] 

event file without randomization 
wavdetect 224.00006(2) -68.72757(1) 

centroid, r = 0.6" 224.00007(2) -68.72757(1) 
ccntroid, r = 2" 224.00002(2) -68.72759(1) 
willi EDEF sul)pi.«>l algorilliiii (>uii)loycd 
ccntroid, r = 0.6" 224.00007(2) -68.72757(1) 
centroid, r = 2" 224.00002(2) -68.72758(1) 



Note. — Count-weighted variances of the ccntroid positions 
are hsted in brackets and refer to the last respective digits. 
These errors do not include astrometric uncertainties. 

angular separation (~ 0".7) from the initial X- 
ray centroid position, we aimed for improving 
the astrometric solution of the event file. As 
reference frames, we used the European South- 
ern Observatory (ESO) New Technology Telescope 
(NTT) SUSI2 observations of PSR B1451-68 (PI: 
Bucciantini, obtained as part of their program 
68.D-0249) and the 2MASS point source catalogue 
(PSC, Skrutskie et al. 2006). The details of the re- 
projection process are given in the Appendix A. 
Our derived overall absolute astrometric uncer- 
tainty of the rcprojccted Chandra event file is 
"^^abs ~ 0.73", and the relative astrometric un- 
certainty of Chandra with respect to the 2MASS 
PSC is Saf'f = 0.66". 

We obtained updated wavdetect positions for 
all X-ray sources on the S3 chip. The X-ray source 
at the pulsar location appears to be point like. 
For this source, we also measured the centroid 
positions of the flux distribution in r = 2" and 
r = 0.6" apertures for both, the event file with just 
the pixel position randomization removed and the 
one with the EDSER subpixel algorithm applied in 
addition. Statistical centroid position errors are 
calculated as the standard error of the mean in 
the pixel coordinates X and Y in the respective 
aperture and transformed to the world coordinate 
system. There are no statistically significant dif- 
ferences between the different positions listed in 
Table 1. 

We extracted a source spectrum within 2" of 



the pulsar ccntroid position. A background spec- 
trum was extracted from an annulus of 5" to 10" 
around the same position. The extracted source 
spectrum counts have been combined to energy 
groups of 15 counts each. Applying XSPEC (ver- 
sion 12.6.0) to analyze the spectrum, we used x^- 
fitting with standard weighting, the tbabs model 
for the X-ray absorption and abundance tables by 
Wilms et al. (2000) as well as photoelectric ab- 
sorption cross sections from Balucinska- Church & 
McCammon (1992) and Yan et al. (1998). Dif- 
ferent spectral models were checked for correspon- 
dence with the energy distribution of the ^ 200 
source counts within an energy range of 0.3 keV 
to 8keV. The fit results are presented in Ta- 
ble 2. We checked whether the VFAINT back- 
ground cleaning influenced the spectral properties 
of the source. Without this correction there are 2 
more source counts (1%), and the corresponding 
spectral model fit results are within la the same 
as listed. 

3. Results and Discussion 

3.1. Position and possible source blending 

In order to determine the pulsar coordinates, 
we combined the Molonglo data used by Sieg- 
man et al. (1993) and Parkes data from 1991 
to 2006. First, we obtained the dispersion mea- 
sure, DM = 8.557 ± 0.014(2<t) cm'^ pc, using the 
Parkes data which have 400, 600 and 1400 MHz 
times of arrivals. This value improves the ac- 
curacy of the previously reported DM = 8.6 ± 
0.2(2(7) cm-3pc by D'Alessandro et al. (1993). 
Fixing the proper motion at the interferomet- 
ric values of Bailes et al. (1990), and using the 
Molonglo as well as the Parkes data, we derive 
the position of PSR B1451-68 as RA(J2000)= 
14:56:00.071(6); Dec(J2000)= -68:43:39.25(5) at 
epoch MJD 50135. The quoted uncertainties 
are twice those given by the least-square solu- 
tion. The expected pulsar coordinates at the 
epoch of the Chandra observations, MJD 55332.7, 
are: RApsr (J2000)= 14:55:59.968(8); DecpsR 
(J2000)= -68:43:39.43(6). The errors reflect 2a 
uncertainties including those due to the proper 
motion uncertainties. 

The measured X-ray source position slightly de- 
pends on the method applied, see Table 1. We will 
use in the following the wavdetect position in the 
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reprojcctccl event file without pixel position ran- 
domization. The measured X-ray source position 
is then RAcxo = 14:56:0.015(7) and Deccxo = 
—68:43:39.25(4). The errors reflect 2cr variances 
of the wavdetect position. 

Close to this position there is a 2MASS 
point source which has also been detected in 
the ESO SUSI2 WB655 image as well as other 
optical surveys. The listed J2000 position of 
2MASS 14560002-6843400 is RAsmass^ 14:56:00.028; 
Dec2MASS= -68:43:40.02 (Cutri et al. 2003). 
2MASS 14560002-6843400 is also listed in the 
USNO Bl catalogue (Monet et al. 2003). Roeser et 
al. (2010) determined the apparent proper motion 
of this source to be —0.6 ± 8.2 mas yr~^ in right 
ascension and —1.9 ± 8.2 mas yr~^ in declination. 
Thus, we can neglect movement of the 2MASS 
source, and the calculated separation between the 
expected pulsar position and the 2MASS point 
source is 0.68" . 

The X-ray source position is offset by 0.26" ± 
0.06" in right ascension and 0.17" ± 0.06" in dec- 
lination from the expected pulsar position. Here, 
the errors reflect only the 3cr wavdetect positional 
variances. In addition, systematic astrometric er- 
rors apply. Our absolute astrometric accuracy is 
better than 3a^^f = 0.73" (Section 2 and Ap- 
pendix A). Thus, the Chandra X-ray position 
is coincident with the expected pulsar position 
within astrometric uncertainties.^ 

The X-ray source position is offset by 
0.07" ± 0.06" in right ascension and 0.76" ± 
0.06" in declination from the catalog position of 
2MASS 14560002-6843400. To discuss the posi- 
tion of this 2A/IASS source, only the Sa^^ll^ = 0.66" 
relative astrometric precision between Chandra 
and 2MASS has to be considered (Section 2 and 
Appendix A). Considering the X-ray statisti- 
cal (wavdetect) positional errors as well as the 
relative astrometric precision between Chandra 
and 2MASS, it appears unlikely that the Chan- 
dra X-ray point source with 0.77" separation 
is the counterpart of the 2MASS point source. 
We discuss spectral constraints on the nature of 
2MASS 14560002-6843400 in the Appendix B.l. 



^The separation between the position of 
2MASS 14560002-6843400 and the optical source de- 
tected in the ESO SUSI2 image is only 46 mas, underlining 
the tight relation between the SUSI2 and 2MASS (and 
subsequently the Chandra) astrometry. 



To estimate a limit of the possible X-ray fhix 
contribution of 2MASS 14560002-6843400 to the 
measured X-ray emission of the pulsar we use dif- 
ferent approaches, which are discussed in detail in 
the Appendix B.2. In short, we measure flux per- 
centages in small apertures centered on the X-ray 
source and the 2MASS source positions, decon- 
volve the image applying the Chandra Ray Tracer 
(ChaRT), MARX and the CIAO task arestore, 
and investigate the count energy distribution of 
the X-ray source. Considering (i) the flux per- 
centage measurements in the reprojected event file 
with or without subpixcl algorithm applied; (ii) 
the comparison of the deconvolved image with an 
deconvolved simulated double source image; (iii) 
the location of the PSF asymmetry region; and 
(iv) the hardness of the detected counts towards 
the 2MASS star; we conclude that any contribu- 
tion of 2MASS 14560002-6843400 to the measured 
X-ray flux of the pulsar must be < 7%. 

3.2. The X-ray spectrum of the pulscir 

As described in Sect. 2, we used XSPEC to 
check different spectral models for the extracted 
X-ray source spectrum. As discussed above, at 
least 93% of the X-ray source flux is attributed to 
the pulsar. We employ in the following the counts 
from the r = 2" extraction region for the pulsar 
spectrum (198 counts in an energy range from 
0.3 keV to 5keV) and neglect the possible (small) 
contribution of the 2MASS star. The fit results are 
presented in Table 2. The derived absorbed X-ray 
fluxes for the different spectral models range from 
2.4 X 10"^^ergcm^s"^ to 2.9 x 10""ergcm^s"^ 
for energies between 0.3 keV and 5keV. Allowing 
all parameters to vary, we tested simple single- 
component models like a BB or a PL spectrum. 
An absorbed, A^h = '2.5tli x 10^^ cm^^ PL model 
with photon index T = 2.4lg;3 flts the spectrum 
well, while the obtained BB fit is statistically unac- 
ceptable. In addition, we checked whether the X- 
ray spectrum can be described by an optically thin 
thermal plasma model as one could expect from a 
stellar corona. The apec plasma model fit is ac- 
ceptable. It is slightly worse than the PL, but bet- 
ter than the BB fit. We discuss the apec param- 
eters with respect to 2MASS 14560002-6843400 
in more detail in Appendix B.l. For the pulsar 
counterpart of the X-ray source, the thin plasma 
model is physically irrelevant and not further dis- 
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Fig. 1.— X-ray spectrum of PSR B1451-68 and 
the best fit absorbed PL model if A^h is a free fit 
parameter; see text and Table 2. The lower panel 
of the upper plot shows the residuals in units of a. 
The lower plot shows 68%, 90%, 99% confidence 
contours in the A'h — T plane of the PL model for 
the pulsar spectral data. 

cussed. 

Among the X-ray spectral fits in Table 2, the 
PL fit and the BB+PL fit show the lowest 
values, but the uncertainties are very large for 
the two-component model. Both fits show sus- 
piciously high best-fit values of the hydrogen col- 
umn density, A^h ~ 2 x fO^^ cm^^. This value is 
close to the Galactic value for this line of sight 
(/ = 313.87°, b = -8.54°): the LAB Survey of 
Galactic neutral hydrogen reports 1.6 x 10^^ cm^^ 
towards this direction (Kalberla et al. 2005), the 
Dickey & Lockman (1990) neutral hydrogen sur- 
vey reports 2.1 x 10^^ cm^^. The pulsar has a 




Energy (keV) 

Fig. 2.— X-ray spectrum of PSR B1451-68 and 
the best-fit absorbed BB-I-PL model in case of 
frozen TVh — 2.6 x 10^°cm~^. The red, solid and 
blue, dashed lines visualize the BB model compo- 
nent and PL model component contributions, re- 
spectively. The lower panel shows the residuals in 
units of a. For more details see text and Table 2. 



parallactic distance of only 450 pc and a low DM 
value. Assuming 10 H atoms for each e~, we de- 
rive an expected = 2.6 x 10^" cm~^ from the 
pulsar DM = 8.557 ± 0.014pccm-3 (Sect. 3.1). 
Freezing A^h at ~ 2.6 x 10^° cm^^, we obtain new 
X-ray spectral fits, presented in Table 2. The PL 
fit is now slightly worse than the BB fit, both are 
unacceptable. The combined BB-I-PL model fits 
the data best. According to the F-test, the com- 
bined model fits the data better than the PL-only 
fit with a probability of 99.4%, and better than the 
blackbody-only fit with a probability of 99.1%. 

The inferred BB emitting area radius is of a 
similar size as in the case of the BB-only fit, 
albeit with larger error. The individual con- 
tributions to the unabsorbed fluxes in the en- 
ergy range from 0.3 keV to 8keV are -Fxwna&s = 



(1.4t°;?i X 

(2.1 ± 1.1) X lO^i'^ergcm^s-i for the BB and 
the PL components, respectively. Flux errors are 
large because the fitted component normalization 
values have high uncertainties due to the low num- 
ber of counts. 

We also checked neutron star (NSA) mod- 
els (Zavlin et al. 1996; Pavlov et al. 1995) with 
fixed parameters for neutron star mass and ra- 
dius (Mns = 1.4Mq, Rns = 12 km). The fits are 



10 -'^"'ergcm^s ^ 



and Fx\,nabs 
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Table 2 
xspec fit results 



Model 


[1020cm-2] 


r 


PL norm at 1 keV 
[10-^photons keV"^ cm-^ s"!] 


kT 
[keV] 


Rbb 
[m] 


EM 
[1038 cm-3] 


red. x^/d-o.f. 


^unabs 
[lO-l^ergcm^s-l] 


PL 

BB 

BB+PL 
APEC 


< 12 

171?? 
6.511;? 


2.2 + 1.0 


13l« 

o 9+10 


'^•^^-o.os 
0.3 + 0.2 

o 1+1.3 


13.2+1122 
Q.Ql^^g^ 


6.6 + 1.0 


0.4/10 
1.4/10 
0.5/8 
0.8/10 


6.0l?i 

2.4+0-^ 

3.8tl:i (PL) /o'7lJ;* (BB) 
0+0.7 

•^•°-0.6 


frozen 


PL 
BB 

BB+PL 


2.6 
2.6 
2.6 


1.7 + 0.2 
1.4 + 0.5 


6.6 + 0.9 

-2.7 


40+0 06 


14.1+1233 

,0 Q+24.2 
ld.»_j2.3 




1.6/11 
1.4/10 
0.5/9 


4.1 + 0.6 

2 4+0 '' 

2.1 + 1.1 (PL) / 1.4l!^;? (BB) 



Note. — All errors indicate 90% confidence intervals, the BB emitting area radius errors taJje the distance error into account. The parallactic distance, corrected 
for the Lutz-Kelker Bias, D = 480l|g pc, is used. The unabsorbed fluxes, F^^nabst given for the energy range from 0.3 keV to 8keV. 



marginally acceptable, and the non-magnetisied 
NSA model fit achieves the lowest x^- However, 
the fits are worse than the BB+PL fit. Allowing 
both the mass and radius to vary results in in- 
plausibly low masses and radii. 

The inferred X-ray spectral fit parameters are 
similar to those found for other old pulsars (e.g., 
Pavlov et al. 2009; Kargaltsev et al. 2006; Becker 
et al. 2004). From a statistical viewpoint alone we 
are unable to differentiate between a purely mag- 
netospheric emission model and a combination of 
magnetospheric and thermal emission. But the 
derived A^h in case of the PL-only fit is a factor 10 
higher than one would expect from the dispersion 
measure. Even accounting for the fit parameter 
error and a possible uncertainty factor of 3 for the 
DM-based A^h, these values do not overlap. There- 
fore, we regard the spectral fits, in particular the 
BB-I-PL fit, with fixed Ah more realistic. In the 
combined model, the photon index of F = 1.4±0.4 
is comparable to the typical values for younger 
pulsars (0.6 < F < 2.1, Kargaltsev & Pavlov 2008; 
Gotthelf 2003), but the inferred BB emission area 
radius, IS.S'ti^'^m is rather small. The conven- 
tional polar cap (PC) radius is Rpc — 280 m as- 
suming a dipolar magnetic field for this pulsar and 
i?NS = 10 km. 

Similar discrepancies have been found for BB 
components of other old pulsars (Pavlov et al. 
2009; Misanovic et al. 2008; Kargaltsev et al. 2006; 
Zhang et al. 2005). Apart from geometrical projec- 
tion effects, the following explanations have been 
discussed: only a small fraction of the PC is heated 
by inflowing particles at the footpoints of "spark 
discharges" created above the PC (Zhang et al. 
2005; Rudcrman & Sutherland 1975); the PCs are 
covered by hydrogen or helium atmosphere rcsiilt- 
ing in effective temperatures of a factor 2 less and 
a radius a factor 3-10 larger than the BB fit (Zavlin 
& Pavlov 2004) ; the nondipolar component of the 
magnetic field is much stronger than the dipolar 
component, causing partially screened accelera- 
tion regions and smaller PCs as consequences (Gil 
et al. 2003). 



3.3. The X-ray efficiency of old rotation- 
powered pulsars 

For PSRB1451-68, the luminosity of the non- 
thermal (PL) X-ray component in the 0.3-8 keV 
energy band is L^**^ = ^Trd^^l^ = 5.9t^;9 x 
10^^ ergs~^, where the error is calculated from the 
90% confidence error of the unabsorbed flux and 
the 80 pc uncertainty of the Lutz-Kelker corrected 
distance. Such a luminosity corresponds to a non- 
thermal X-ray efficiency ??nonth = i^"'''^/-^ ~ 
3 X 10~^. The X-ray luminosities versus character- 
istic ages of old pulsars, including PSRB1451-68, 
are plotted in Figinc 3. There is no obvious char- 
acteristic age effect visible in this plot. We also 
include PSR B1451-68 in the plot showing the X- 
ray luminosities versus spin-down energies of old 
pulsars (Figure 4). Distances to other old pul- 
sars are updated according to Deller et al. (2009) 
and corrected for the Lutz-Kelker bias (Verbiest 
et al. 2010). The non-thermal X-ray efficiency of 
PSR B1451— 68 is larger than those of yoimg pul- 
sars, most of which have r/nonth values smaller than 
10-3, see Figure 5. The efficiency of PSRB1451- 
68 is, however, comparable to those of other old 
pulsars (Fig. 4). While it is possible that old pul- 
sar spectra with too few counts may have an un- 
detected thermal contribution (e.g., from heated 
polar caps). Figure 4 shows only little influence 
of such a thermal component on the inferred non- 
thermal X-ray efficiencies in the energy range 1 
to lOkeV for the 3 pulsars with sufficiently large 
number of counts. This supports the hypothesis 
that the conversion of spin-down power into X-ray 
emission becomes more efficient as pulsars get less 
powerful (Zharikov et al. 2006; Kargaltsev et al. 
2006). There is, however, an observational bias in 
favor of brighter X-ray sources which can lead to 
the detection of only the most efficient old neutron 
stars. 

A look on the increasingly complete observa- 
tional statistics of the closest pulsars can help 
to minimize the effect of the observational bias 
in favor of brighter X-ray sources. Instead of 
the characteristic ages, which may be very dif- 
ferent from the true ages, we prefer to use the 
directly measurable E as parameter to differen- 
tiate between powerful, young and less powerful, 
old pulsars. We choose E = 9 x lO^^ergs"^ as 
a rough, admittedly arbitrary boundary between 
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young and old pulsars. This boundary was mo- 
tivated by pulsar statistics in the ATNF pulsar 
catalog^ (Manchester et al. 2005), where most ra- 
dio pulsars have characteristic ages t > 1 Myr for 
E < 9 X 10'^^ erg s"'^ while younger radio pulsars 
usually have higher E. The close 7-ray pulsar 
PSR J1741-2054 with r = 386 kyr, but modest 
E = 9.5x 10'^^ ergs~^ influenced our choice of the 
particular boundary value. 

According to the ATNF pulsar catalog, PSR Bf 45f- 
68 with its Lutz-Kelker corrected distance of 
4801 gopc is the eleventh closest of the isolated 
radio pulsars with ii^ < 9 x lO'^'^ ergs^^. Six 
of these eleven pulsars'^ are now detected in X- 
rays, all six have parallactic distances. Two of 
the X-ray detected pulsars, PSRB1929-flO and 
PSRB0823+26, have X-ray efficiencies lO"-* < 
?7nonth < 10^'^, the other four X-ray detected pul- 
sars have 77nonth > 10^'^. Thus, at least 36% of all 
the closest old pulsars with E <9 x 10"^'^ ergss"^ 

show r?nonth > 10"^. 

For comparison, we consider now the 11 clos- 
est, isolated, younger pulsars with £■ > 9 x 
10'^^ ergs s~^. Nine of them are detected in X- 
rays, the two others were not probed deep enough. 
Four of these nine detected pulsars have paral- 
lactic distances. Based on information of the 
ATNF pulsar catalog. X-ray flux investigations by 
Marelh et al. (2011); Camilo et al. (2009); Kargalt- 
sev & Pavlov (2008); Tepedelenlioglu & Ogelman 
(2007), and distance corrections according to Ver- 
biest et al. (2010) and Mignani et al. (2010), we 
find that all nine X-ray detected younger pulsars 
have T^nonth < 10^"^. Thus, at most 18% of all 
these 11 pulsars with _E > 9 x lO'^'^ ergs s~^ could 

have 77nonth > 10"^. 

Thus, looking only at the 11 closest sources 
in each case, it is suggestive, that at least one 
third of old pulsars have efficiencies rj > 0.001, 
while none of the detected closest, younger pul- 
sars has. Recently, Marelli et al. (2011) reported 
the following Lx-E relation for 29 Fermi pulsars 
with available distance estimates: logj^g Lx,29 = 



(1.11 



+0.21-. 
0.30^ 



(1.04 ± 0.09)logio£^34, which is in 



agreement with older estimates for X-ray pulsars 



^ http: / / www. atnf . csiro.au / research /pulsar /psrcat/ 
^PSR B2224-I-65 with its potentially overestimated X-ray 
luminosity due to an unresolved shocked pulsar wind com- 
ponent is not among the eleven closest old pulsars. 



by Kargaltsev & Pavlov (2008) and Possenti et 
al. (2002). From this, one would indeed expect 
ry < 10"'^ for most X-ray pulsars, and the > 36% 
deviation of the closest, low E pulsars appear 
even more puzzling. Unfortunately, an insufficient 
number of probed objects, inhomogeneous depth 
of the X-ray observations, low signal-to-noise ra- 
tios in many cases, potentially unresolved compact 
pulsar wind nebulae, as well as large uncertainties 
of DM-based distances, prohibit currently any fur- 
ther conclusive statistics, for a larger sample of old 
pulsars in particular. Whether there is indeed a 
higher X-ray efficiency for pulsars with low E can 
only be checked by increasing the sample of such 
pulsars with sufficiently sensitive X-ray observa- 
tions. 
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Fig. 3. — X-ray luminosities and upper limits of 
eleven old pulsars versus their characteristic ages. 
In case of a sufficient number of counts, the di- 
amonds and asterisks show the non-thermal and 
thermal X-ray luminosities, respectively; other- 
wise the luminosities were obtained from PL fits. 
The arrows mark upper limits derived from X-ray 
non-detections. This figure is an update of the one 
presented earlier by Kargaltsev et al. (2006) using 
new distances by Deller et al. (2009), the Lutz- 
Kelker bias corrections by Verbiest et al. (2010), 
and the upper hmit resuh for PSR J2144-3933 
by Tiengo et al. (2011). Note that the large X- 
ray luminosity of the extremely fast moving PSR 
B2224-I-65 may be overestimated due to a sus- 
pected, unresolved shocked pulsar wind compo- 
nent (e.g., Johnson & Wang 2010). 
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Fig. 4. — X-ray luminosities and upper limits of 
eleven old pulsars versus their spin-down power. 
In case of a sufficient number of counts, the di- 
amonds and asterisks show the non-thermal and 
thermal X-ray luminosities, respectively; other- 
wise the luminosities were obtained from PL fits. 
The arrows mark upper limits derived from X- 
ray non-detections. This figure is an update of 
those presented earlier by Kargaltsev et al. (2006) 
and Pavlov et al. (2009) using new distances by 
Deller et al. (2009), the Lutz-Kelker bias correc- 
tions by Verbiest et al. (2010), and the upper limit 
result for PSR J2144-3933 by Tiengo et al. (2011). 
Note that the large X-ray luminosity of the ex- 
tremely fast moving PSR B2224-I-65 may be over- 
estimated due to a suspected, unresolved shocked 
pulsar wind component (e.g., Johnson & Wang 
2010). 



factor 20 less than expected conventionally. The 
nonthermal X-ray efficiency, 3 x 10~^, is high 
in comparison to younger pulsars, but comparable 
to estimates for other old pulsars. 
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4. Summary 

We investigated Chandra observations of PSR 
B1451— 68. Through various image analysis tech- 
niques we estimate the possible contribution from 
a nearby 2MASS star to the pulsar's X-ray flux 
to be less than 7%. No significant extended emis- 
sion is seen. The pulsar has a soft X-ray spec- 
trum similar to those of other old pulsars. Nom- 
inally, the spectrum is best fit with a power law, 
having a photon index of F = 2.41q'3. How- 
ever, the inferred hydrogen column density, iVn, 
is uncomfortably close to the Galactic value of HI 
in this direction. Fixing iVn to the DM-derived 
2.6 X 10^°cm~^, a combination of a thermal and 
a non-thermal component fits the data best. The 
inferred radius of a heated polar cap is small, a 



9 



36 - 



> 

(D 

O 
I 



o 



34 



32 - 



30 -X 



28 



❖ - 



- i> 



I- 




32 



34 



36 



38 



Log E 



Fig. 5. — Non-thermal X-ray pulsar luminosities and upper limits versus the pulsar spin-down power. This 
plot is an update of Figure 5 by Kargaltsev & Pavlov (2008) including the old pulsars from the previous 
figure, recycled X-ray pulsars and new X-ray detected Fermi pulsars. The latter will be presented in detail 
in Kargaltsev et al. 2012 (submitted). As in the previous figure, parallactic distances were corrected for the 
Lutz-Kelker bias, where applicable. 
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A. Reprojection of the X-ray data and X-ray source position 

For the optical reference frame, we used the European Southern Observatory (ESO) New Technology 
Telescope (NTT) observations of PSR B1451-68. The NTT was equipped with the SUSI2 imager. We chose 
an image which was obtained with the WB655, a wide ii-band filter in February 2002. The field of view 
of the CCD chip with the target on it is 2.8' x 5.5', and there are plenty of 2MASS point sources in this 
field, e.g., 46 with the highest quality flag AAA. We obtain the absolute astrometry of the ESO SUSI2 image 
with the help of the 2MASS point source catalog (PSC, Skrutskie et al. 2006) using Graphical Astronomy 
and Image Analysis Tool (GAIA) (Draper et al. 2007). Nominally, an absolute positional accuracy of 
3cr2MASS = 300 mas (Skrutskie et al. 2006) can be achieved for high-quality 2MASS point sources. The rms 
of the GAIA astrometric fit to the ^^^-2MASS source positions was determined to be <tsusi2 = 178 mas. 
Thus, we expect the source positions in the ESO SUSI2 image to have an absolute astrometric 3a error of 
612 mas. 

Nine optical sources in the SUSI2 image were found to correspond to wavdetect X-ray sources on the ACIS 
S3 chip. Five of them are also 2MASS point sources. We excluded one source close to the expected pulsar 
position and proceeded with eight optical detections. The sources are indicated in Figure 6. The X-ray 
and optical sources were matched and the Chandra image astrometry updated by applying the CIAO tasks 
reproject_aspect and wcsjnatch. After removing one poor match, seven sources remained. For the seven 
sources the average residual after reprojecting was acxo = 130 mas in comparison to 340 mas for the original 
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Fig. 6. Left: the NTT SUSI2 WB655 image is shown for the CCD chip covering the target PSR B1451-68. 
Right: the reprojected X-ray image is plotted for the same coordinate range. The approximate size of each 
image is 2.5' x 3.8'. The brightness scale is streched to highlight the respective sources. The X-ray source 
position of PSR B1451— 68 is marked by the largest box in the upper part of the images. Smaller boxes in 
the SUSI2 image indicate optical sources used for the reprojection. The circles in the X-ray image indicate 
the corresponding X-ray sources. 
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event file (all without pixel position radomization) . Thus, our overall absolute astromctric uncertainty is 
3 (^iusi2 + '''Imass + ^cxo)^ = 3(T^J^ = 0.73", while the relative astrometric uncertainty of Chandra with 
respect to the 2MASS PSC is 3 (ct^usk + c^cxo)^ = ^'^ref = 0-66"- 

B. 2MASS 14560002-6843400 

B.l. Spectral constraints on the nature of 2MASS 14560002-6843400 

Here, we want to test if neglecting positional arguments - 2MASS 14560002-6843400 could be in prin- 
ciple the counterpart of the X-ray source. As described in Sect. 2, we checked different spectral models for 
correspondence with the energy distribution of the ^ 200 source counts within an energy range of 0.3 keV to 
8keV; the fit results are listed in Table 2. Amongst others, we checked whether the X-ray spectrum can be 
described by an optically thin thermal plasma model as one could expect from a stellar corona. In Table 2, 
we list the obtained apec model parameter values, the mekal or raymond models gave similar values. The 
apec plasma model fit is acceptable. The derived plasma temperature is high, kT = 3.1^Q'gkeV. Such high 
temperatures can be reached in young stars (e.g., Preibisch et al. 2005; Getman et al. 2005). 
As a reminder, the derived absorbed X-ray fluxes for the different spectral models range from 2.4 x 
10~^^ergcm^ s~^ to 2.9 x 10~^^ergcm^ s~^ for energies between 0.3keV and 5keV. We use the formula 
by Maccacaro et al. (1988) log(/x//y) = log(/x) -I- 0.4TOy -|- 5.37 to estimate the optical to X-ray flux ratio. 
2MASS 14560002-6843400 has an apparent l/-band magnitude of my = 16.2 mag in the General Guide Star 
Catalogue version 2.3.2 (Lasker et al. 2008), and my = 15.7 mag in the NOMAD catalogue (Zacharias et 
al. 2004). Assuming that the detected X-ray source is the counterpart of 2MASS 14560002-6843400, the 
log(/x//y) range is accordingly from —1.7 to —2.0. Such flux ratios are typical for K or M stars (e.g., 
Agiieros et al. 2009). 

2MASS 14560002-6843400 has the following magnitude measurements: J= 14.50 ± 0.04 mag, H= 13.94 ± 
0.04mag, K= 13.83 ± 0.06 mag, B2= 16.83mag (USNO Bl), R2= 15.05 mag (USNO Bl), /= 15.11mag 
(USNO Bl). The obtained colors B2 — V and V — I, are consistent with a late G to early K star, the colors 
J — H, H — K, V — J, V — H, and V — K indicate a KO to K2 star following intrinsic stellar colors listed by 
Currie et al. (2010) and Covey et al. (2007). A K2 dwarf would be consistent with all colors without requiring 
significant extinction. From the Nn = 6.5^47 x lO^'' cm^'^ obtained for the apec model fit we estimate the 
optical extinction Ay = O.Sd'^Q ^e using the relation between the hydrogen column density and extinction 
by Predchl & Schmitt (1995). Using the extincion relations by Cardelli et al. (1989) and Ay = 0.36 wc 
estimate color corrections and find that a G9 dwarf is the most likely extinction-corrected counterpart for 
the cataloged NIR and optical magnitudes. Assuming Mj = 4.2 mag according to the Padova tracks (Bertelli 
et al. 2008), the reported apparent J magnitude translates into a distance of 1.1 kpc. For a K2 dwarf without 
extinction, assuming Mj = 4.8 mag, the reported apparent J magnitude translates into a distance of 0.9 kpc. 

A typical K or late G dwarf main sequence star is expected to have its X-ray peak emission below 1 keV 
and show nearly no emission above 1 keV. Young dwarf stars, on the other hand, could produce the observed 
X-ray spectrum. There is no obvious star formation region within 1 degree of the target position, and as 
noted above, 2MASS 14560002-6843400 has a very slow proper motion (< 12masyr~^). While one cannot 
exclude an in-situ young stellar object, there is no indication for this from extinction, outflows, or flares, and 
a late G or early K dwarf main sequence star seems to be the most likely counterpart of 2MASS 14560002- 
6843400. Prom this in turn, we conclude that the contribution of 2MASS 14560002-6843400 to the X-ray 
emission is small, and that the X-ray source is not the counterpart of 2MASS 14560002-6843400, but of the 
pulsar. 
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Fig. 7. — This RGB image shows the reprojected event file with the CIAO EDSER subpixel algorithm applied. 
A 0.25 subpixel binning, as well as slight smoothing (Gauss kernel 2 subpixel), was used in each band. The 
red band shows the events from 0.3 keV to IkeV, the green band shows the events from IkeV to 2keV, the 
blue band shows the events from 2keV to 8keV. The scale is linear and the same for all bands. North is up, 
East is to the left. The expected pulsar position is marked with a black cross, the 2MASS catalog position 
of a nearby star is marked by a red diamond, the position of the star in the SUSI2 image is marked with 
the red square. The big cyan circle has a radius of 0.6" and marks the centroid position of the X-ray source. 
The other cyan region indicates the area of the Chandra point spread function asymmetry as inferred by 
applying the CIAO task make_psf _asjnninetry_region. See text for the discussion. 
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B.2. Investigation of possible source blending 

Wc wish to obtain an \ippcr limit for any potential X-ray flux contribution of the 2MASS star to the 
X-ray source detected at the pulsar position. First, we look at the spatial distribution of the dominant 
count energies in the X-ray image. Figure B.2 shows a color-coded image of the reprojected event file 
with the CIAO EDSER subpixcl algorithm (Li et al. 2004) applied. Apparently, there is emission near the 
2MASS source. However, the location of the 2MASS soure is very close to a region where enhanced emission 
is expected due to the recently found asymmetry in the Chandra point spread function^ (PSF). For our 
observation wc obtained the inflicted region by applying the CIAO task make_psf .asymmetry jregion, it is 
marked in the figure in cyan and the artifact can constitute ~ 5% of the total brightness. If a normal K 
star at the position of the 2MASS object contributes to the X-ray flux, one would expect an enhancement 
in soft counts towards this position. However, there appear to be slightly more hard counts in the direction 
of the 2MASS source. A hard spectrum is unusual for a normal star, less so for a young stellar object. As 
discussed in Appendix B.l, the 2MASS source is unlikely to be a young stellar object. 

Next, we measure the number of counts in similar regions around the center of the X-ray source to quantify 
the overabundance of counts towards the 2MASS star. Figure 8 shows the chosen counting regions: The 'star' 
circle is centered on the 2MASS source's 2002 optical SUSI2 position. Circles 1 to 4 have similar separations 
from the X-ray centroid and the same size as the 'star' region. They are used to obtain the average expected 
number of counts from the PSF wings for this particular separation and area. This assumes a symmetric 
PSF (see comments below). To relate the 'star' counts, Cstar- to the pulsar counts, CpsR. we use a circle 
of the same size at the position of the X-ray centroid. We estimate the ratio between the extra 'star' 
counts and the pulsar counts as (Cstar — (Ci + C2 + C3 + C4) 0.25)/Cpsr. Of course, this is a very crude 
estimation neglecting a detailcid Chandra PSF model in general and the outer PSF wings in particular. In the 
ideal case of well separated point sources with similar spectra and using the same aperture for each source, 
their flux ratios should be constant for different aperture sizes. In our case we need to account for a flux 
contribution due to the overlapping of the individual PSFs. This contribution is significant for the potcintial 
'star' source, hence the subtraction of 0.25 (Ci + C2 + Cs + C4) is an approximation of this contribution. 
The flux contribution of the 'star' PSF to Cpsr is negligible as seen from Chandra MARX simulations. For 
circles with r = 0.4" wc infer 'star' percentages of 6.3% for the event file with the EDSER subpixel algorithm 
applied, and 0% (actually -1.3%) for the event flle without randomization. For circles with r = 0.3" we 
infer 'star' percentages of 7.1% for the event flle with the EDSER subpixel algorithm applied, and 3.1% for 
the event file without randomization. Note that the ACIS encircled energy is 50% for an aperture radius 
of r = 0.418", while it is 90% for r = 2" , the radius used in Sections 2 and 3.2 {Chandra Proposers' 
Observatory Guide, sect. 6.6 ^). Due to the Poisson noise character of the very few counts we are dealing 
with, it is difficult to obtain a rigorous error estimate for the above formula. Changing slightly the positions 
of the individual counting regions, we approximate the error to be around 1% of the pulsar counts. The 
known PSF assymmetry region of enhanced flux is actually partly overlapping our 'star' circle (see Fig. 8). 
Since part of the enhanced flux in the 'star' circle is likely due to the asymmetry 'leakage', our estimated 
percentages are indeed conservative upper limits. 

We also deconvolved the image applying the Chandra Ray Tracer (ChaRT), MARX (v. 4.4 and 4.5) 
and the CIAO task arestore following standard CIAO threads. Since the subpixel algorithm is not yet 
implemented in the simulators, we discuss here only the event file without pixel radomization. We apply 
ChaRT (Carter et al. 2003) to construct the PSF at the source position. As the input spectrum, we used 
the parameters derived from the PL fit of the source spectrum (see Sect. 3.2) in an energy range from 
0.3 to SkeV. From the centroid source position we obtain the off-axis angle, 9 = 18". 6, and the mirror 
spherical coordinate azimuth, — 327°. 9^. We use a long exposure time (180 ks) in order to generate a high 



^http: / /cxc. harvard.edu/ciao/caveats/psf_artifact. html 
^http://cxc. harvard.edu/proposer/POG/html/index.html 

^For details on the Chandra coordinate systrems see http://cxc.haj-vaj-d.edu/contrib/jcm/ncoords.ps 
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signal-to-noisc-ratio PSF image which is advantageous for the following deconvolution and the souree extent 
estimate. The raytrace file is then used within MARX*" to produce a PSF image with 0.25 pixel binning. 
We set the DitherBlur parameter to 0.2" as recommended in the MARX manual for ACIS data without 
pixel randomisation. Using the PSF image, the CIAO task srcextent reports the X-ray source to be not 
extended at 90% confidence. We then deconvolved the event file applying the task arestore which is based 
on the Lucy- Richardson deconvolution algorithm (Lucy 1974). The image is shown in Figure 9. The new 
centroid position of the main source in the deconvolved image has a small offset of 0.1" with respect to 
the centroid position in the original image. Comparing the count numbers in a r = 0.4" circle centered 
on the former and a circle of the same size centered on the SUSI2 optical position of the 2MASS star, we 
derive a flux percentage of 6.5% for the counts around the star with respect to those of the pulsar. Using 
the original centroid position for the main source decreases this percentage to 6.1%. The deconvolution 
produces 'extra' counts not only towards the direction of the 2MASS source, which is located at one edge 
of the PSF asymmetry region, but also southwest of the pulsar where the other edge of the PSF asymmetry 
region lies. Since neither ChaRT, MARX, nor arestore account for this asymmetry, some artefacts in the 
south-southwest direction can be expected after deconvolution. 

For comparison, we also used MARX to simulate two sources having combined as many counts as the r = 2" 
main source region - one source with 95% of the flux at the main source centroid position, and one with 5% 
flux at the SUSI2 optical position of the 2MASS star. Again, we set the DitherBlur parameter to 0.2" and 
proceeded with the processing as described above. Again we find an ^ 0.1" shift (but in another direction) 
in the centroid position of the deconvolved image comparing to the original input coordinates. While these 
shifts are all within the nominal Chandra astrometric accuracy, we do not utilize the arestore positions for 
further astrometric purposes. Using again the r = 0.4" circles as indicated in Figure 9, we find 4.0% of the 
simulated main source counts in the circle around the simulated 2MASS source. From this we infer that, in 
principle, it is possible to recover a ~ 5% X-ray source at a 0.77" separation, which the 2MASS source has 
from the main X-ray source centroid. 



^http://space/mit/edu/CXC/MARX/index.html 
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Fig. 8. — This image shows same-size source regions on the reprojected event file with the CIAO EDSER 
subpixel algorithm applied. The binning is 0.25 Chandra ACIS sky pixel. The 'star' circle is centered on 
the SUSI2 optical position of the 2MASS point source - its center position is marked by the black box. 
Similarly, the white box marks the centroid position of the X-ray source. All circles have radii of 0.4". 
In the southwest, the region of the Chandra PSF asymmetry is indicated as obtained with the CIAO task 
make_psf -asjnninetryjregion. See text for the discussion. 
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Fig. 9. — On the left is the deconvolved image produced from the event file without randomization, obtained 
by applying the CIAO task arestore. The right image is the deconvolved image, which was produced from 
a simulated event file with 2 sources - one with 95% flux and the other with 5% flux of the detected X-ray 
source in our data set. The dashed circle with r = 0.6" is centered at the centroid position of the original 
event file of our observation. The pulsar position is marked with a cross, the 2MASS point source position 
with a diamond, the 2002 SUSI2 optical position of this 2MASS source with a small box. The region of the 
PSF asymmetry is also indicated. One of the two solid, r — 0.4", circles in each image is centered on the 
corresponding centroid position of the deconvolved main source, the other is centered on the SUSI2 position. 
Both images are in log brightness scale to emphasize the low count numbers. See text for the discussion. 
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